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X ■ Abstract 

* ^ W reeen, th^oop reSU Us „» th e qUark a„ d glu o„ ^ tu „e tl o„s a„ d form f a Ct0rS , 
we derive the complete threshold-enhanced third-order (N 3 LO) QCD corrections to the total cross 
sections for the production of lepton pairs and the Higgs boson in hadron collisions. These results, 
for the latter case obtained in the heavy top-quark limit, are employed to extend the threshold 
resummation for these processes to the fourth logarithmic order. We investigate the numerical 
impact of the higher-order corrections for Higgs boson production at the Tevatron and the LHC. 
Our results, suitably treated in Mellin A^-space, provide a sufficiently accurate approximation to 
the full N 3 LO contributions. Corrections of about 5% at the LHC and 10% at the TEVATRON are 
found for typical Higgs masses. The N 3 LO predictions exhibit a considerably reduced dependence 
on the renormalization scale with, for the first time, stationary points close to the Higgs mass. 



The production of lepton pairs and especially the Higgs boson H, seeRefs. [1] for detailed reviews, 
are among the most important processes in high-energy proton collisions. The corresponding 
cross sections receive sizable higher-order QCD corrections, necessitating calculations beyond the 
standard next-to-leading order (NLO) of perturbative QCD. After the early calculation of the next- 
to-next-to-leading order (NNLO) corrections to the Drell-Yan process in Refs. [2,3], considerable 
progress has been achieved in this field during the last five years. The NNLO corrections are now 
completely known also for Higgs production in the heavy top-quark limit [4-10], and the all-order 
resummation of the threshold logarithms [11,12] has been extended to the next-to-next-to-leading 
logarithmic (NNLL) accuracy for both processes [13, 14]. Still, at colliders energies especially for 
Higgs production, corrections of yet higher order are not entirely negligible. 

In this letter we derive the complete logarithmically enhanced soft-emission corrections to both 
lepton pair and Higgs boson production at the third order (N 3 LO) in the strong coupling constant 
cc s , and extend the corresponding threshold resummations to the N 3 LL contributions. The first 
step is achieved by analysing the general mass-factorization structure of the third-order coefficient 
functions c 3 in terms of recent results on the three-loop splitting functions [15, 16] and the quark 
and gluon form factors [17, 18], analogous to the second-order procedure of Ref. [19]. For the 
second step we just have to determine one resummation coefficient, usually denoted by D$, from 
these results, as the structure of the resummation exponent has already been derived in Ref. [20] to 
the required accuracy. As discussed below our results on c 3 for Higgs production, suitably treated 
in Mellin /V-space, provide a very good approximation to the complete N 3 LO corrections, thus 
facilitating improved predictions for the cross sections at the TEVATRON and the LHC. 

In the soft limit, i.e., retaining only contributions of the forms 
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to the coefficient functions, only the respective subprocesses qq — > y* — > l + l~ and gg — > H 
contribute to the Drell-Yan process and Higgs boson production. In the latter case, the Hgg vertex 
is an effective interaction in the limit of a heavy top quark, 

£eff= —ChHG^G"^ , (2) 

where G^ v denotes the gluon field strength tensor, and the prefactor Ch includes all QCD correc- 
tions to the top-quark loop. This coefficient is of order cc s and known up to N 3 LO (oc^) [21]. The 
analysis of the higher-order corrections in the heavy-top limit is justified by the agreement between 
this approximation and the full calculations at NLO [22-24]. 

The general structure of the expansion coefficients W„ of the bare partonic cross section, 



n=0 

is given by 



£ (affw^ a s =|, (3) 
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W$ = 8(1 -x) 

Wf = 2ReJiS(l-x)+5i 

W 2 b = (2Re^ 2 + |Ji| 2 )5(l-^) + 2Rej 1 5 1 +5 2 

W 3 b = (2Re^ 3 + 2|^iiF2|)5(l-x) + (2ReiF2 + |Ji| 2 )5i+2ReJi52 + 53. (4) 

Here J„ denotes the bare n-loop time-like quark or gluon form factor, calculated in dimensional 
regularization with D = 4 — 2e and, as all quantities, expanded in terms of a s = oc s / (4ti). The 
dependence of the pure real-emission contributions S n on the scaling variable x = M 2 /s is given 
by the D-dimensional +-distributions fin e defined by 

f k 4x) = e[(l-x)- l - k£ ] + = -l5(l-*) + £^^£^ (5) 

with of Eq. ©. As appropriate for a parallel tretament of the two processes, the coefficient 
function for Higgs production is defined as in Ref. [8], i.e., Ch in Eq. © is kept as a prefactor. 

The expansion coefficients Wk obtained from Eq. (0) after renormalizing the coupling constant, 
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a b = as _^ + , E£_^± + 5 (6) 



and multiplying, for Higgs production, with the square of the renormalization constant [25, 26] 

Z G 2 = [l-p(a s )/(a s £)]- 1 , P(ocs) = -(3 a s 2 -(3 1 a s 3 -(3 2 a s 4 -... , (7) 

of the operator G^ v G a '^ v in Eq. © obey the Kinoshita-Lee-Nauenberg theorem [27,28] and the 
mass-factorization relations (putting fj 2 = M 2 * H ) 

Wi = -Yo + cx+eaj + ^j, (8) 
W 2 = ^{(2Yo-Po)Yo} + -{Yi+2c 1 Yo} + c 2 + 2fl 1 Yo + e{«2 + 2fe 1 Yo} ) (9) 

W3 = ^3-{2(Yo-Po)(2Yo-Po)Yo} + ^2{6YiYo-2poYi-2piYo 

+ 3cj (2 Y - Po) Yo) + ^ {2Y 2 + 3q Yi + 6c 2 Yo + 3a, (2 Yo - (3 ) Yo) 

+ C3+a 1 Yi+2a 2 Yo + ^i(2Yo-Po)Yo- (10) 

Here the anomalous dimensions Y^ are related by a conventional sign to the diagonal Altarelli-Parisi 
splitting functions P^. In x- space, where these quantities (in the MS scheme adopted throughout 
this article) have soft limits of the form [29] 

P k , = A k V +P k & 8(l-x) , (11) 
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the products in Eq. ([8]) - (1101) have to be read as Mellin convolutions. To the required accuracy 
these convolutions can be readily carried out using, for example, the appendix of Ref. [30]. 



Eqs. dH]) - (|T0l) can be used to derive all contributions to the MS coefficient function c n , 
once the coefficients A n , P„ are known together with all 1/e pole terms of the n-loop form factor 
and suitable lower-order information. The salient point for this extraction is the structure © of 
the soft emissions linking the coefficients of £ 8(1 — x) to those of (Dq\ thus a mass-factorization 
constraint on the former term fixes the latter coefficient. With the results of Refs. [2, 8, 9] and 
[15-18] the above conditions are fulfilled at n = 3 for both lepton pair and Higgs boson production. 

Thus, treating Higgs boson production first, we insert the £-expansion 



S n = fin.e £ 2ns nl e l (12) 

l=-2n 

into Eqs. ([8]) - (flOl) and recursively determine the coefficients s The first-order result is known 
to all orders in £. For later convenience, we here present its expansion up to order 8 3 , 
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The corresponding second-order coefficients read 
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Here denotes the number of effectively massless quark flavours, Cp and C A are the usual SU(N) 
colour factors, with Cp = 4/3 and Ca = 3 for QCD, and C, n represents Riemann's ^-function. 

From 5 2 _ 4 . . . s 2 -i an( l me first-order results one recovers the terms [4, 5] of the NNLO 



coefficient function at the renormalization and factorization scales /j 
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The SU(N) result (fT4|) for the coefficient s 2 q, on the other hand, is fixed by the corresponding 
5(1 — x) contribution to Eq. (fT51) derived in Ref. [9]. Actually this colour-factor decomposition is 
checked (and could have been predicted from the N = 3 QCD results of Refs. [4,5]) by the absence 
of a CpUj term in S2 obvious from the colour structure of the contributing Feynman diagrams. 

The quantity a 2 in Eq. ® has not been computed so far, hence the coefficient s 2 1 is unknown 
at this point. This suggests a problem, as this coefficient enters the £° part of Eq. (flOl) . However, its 
contribution to the CDq term (but not the 8(1 —x) piece) of c 3 is found to cancel in the end. Keeping 
s 2 l as an unknown in the intermediate relations, the third-order coefficients in Eq. (fT3l are 
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Analogous to s 2 discussed above, the coefficient s 3 cannot be derived by mass -factorization 
arguments, but requires a third-order calculation like in the case of deep-inelastic scattering [31]. 

The above results, after combination with the gluon splitting function [16] and form factor [18] 
according to Eqs. © and ©, lead to the following soft-emission contribution to the third-order 
(N 3 LO) coefficient function for Higgs boson production at /j 2 = /j 2 = M^: 
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Eqs. (fTTT) - (12TI) agree with the results derived from the NNLL threshold resummation in Ref. [14]. 
Eq. (|22l) represents a new result of the present study. 

We now turn to lepton pair production. The corresponding coefficients s n l can be derived 
in the same manner, some being related to their Higgs counterparts by simple replacements of 
Cf by Ca- Specifically, S\ for the Drell-Yan case is obtained from Eq. (fT3l by a factor Cf/Ca, 
and s n _ 2n+ k f° r n>2 and k = 0, 1 by multiplication with (Cf /Ca)* from Eqs. dT4|) and (fT6l) . 
Accordingly the coefficients of ©5 and ©4 of the third-order coefficient function for the Drell- 
Yan process are related to Eqs. (fTTT) and (fT8~1) . respectively, by factors C F /C A and C F /C A . The 
remaining coefficients for the standard scale choice // 2 = fij = M 2 * are given by 
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Also Eqs. (1231) - (1231) agree with the results derived from the NNLL threshold resummation [13], 
while Eq. (1261) is a new result of this study. The additional coefficients for /jj ^ A/ 2 * and // 2 ^ 
nj can be derived analogously to Eqs. (2.16) - (2.18) of Ref. [32] or using the N 3 LL threshold 
resummation expression [20], but will be skipped here for brevity. 

We note that the ^-function terms of highest transcendentality n, i.e, the coefficients of C, n and 
£i £j with i + j = n, in the £~ 2 '+" contributions to the pure real-emission function Si agree between 
Higgs production and the Drell-Yan process for the Super- Yang-Mills case C A = C F = n c . The 
same holds for the quark and gluon form factors [17, 18] and, consequently, also for the ^-function 
terms of weight n in the soft logarithms 2) 2 z-i-n °f me coefficient functions for Higgs boson and 
lepton pair production, see Eqs. (fT71)- (|22l) and (T23l-(l26l). By construction, generalizing Eq. ©, 
this feature extends to all orders of perturbation theory. 

We now turn to the threshold resummation. For the processes under consideration, the coeffi- 
cient functions exponentiate after transformation to Mellin /V-space [11,12], 
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Here g 0k collects the N- independent contributions at the k-th order, and the resummation exponent 
G N takes the form 



J N (Q 2 ) = \nN- gl (l)+g 2 (l)+a s g 3 (l)+a 2 g 4 (l) 
'. The functions g-, and g A have been determined in 
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with X = |3oa s In N. The functions g 3 and g 4 have been determined in Refs. [13, 14] and [20], 
respectively. Besides the quantities in Eq. (fTTI) and lower-order coefficients, the functions gk 
depend on one parameter, usually denoted by \. 

Before we present our new results for the coefficient D3, we recall, for the convenience of the 
reader, the /V-independent first- and second-order contributions which enter its determination. For 
Higgs boson production these coefficients read 



g 0l = C A (16^ 2 + 8y e 2 ) 
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The corresponding results for the Drell-Yan case are given by 
g i = C F (-16 + 16^ + 8y, 
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Inserting the above results into the explicit formulae for the resummation exponents [13, 14,20], 
we recover the known coefficients (as always referring to the expansion parameter a s = oc s / (4%) ) 
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and derive the new third-order contribution 
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with Cj = Cf for the Drell-Yan case, and Cj = C A for Higgs production. Hence, not unexpectedly, 
we find that also D3 is maximally non-abelian, with the quark and gluon cases related by an overall 
factor Cp/C A . This is the same behaviour as shown by the cusp anomalous dimensions A n in 
Eq. (fTTI) [29] and by the form-factor resummation coefficients f n known up to three loops [9, 18]. In 
fact, there is a simple relation between the coefficients D n and /„ (using the notation of Ref. [18]), 



D 2 = 2pVi, -2/2 

D 3 = 2p 1 5 1 , -4p^ M +4po(52,o- 3 6/5C2 2 C / 2 )-2/ 3 



( 3 6) 



of which the first line of has been derived before in Ref. [33] from an extension of the threshold 
resummation to /^-independent contributions. In the present mass-factorization framework, the s n l 
terms in Eqs. (l36l) can be traced back to the coupling-constant renormalization of Eqs. ©. 
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We now turn to the numerical impact of the N LO and resummation corrections to the coef- 
ficient functions, confining ourselves to Higgs boson production for brevity. All parameters are 
taken over from Ref. [8], i.e., we use the values m t = 173.4 GeV (very close to the present world 
average) and Gf = 4541.7 pb for the top mass and Fermi constant in the prefactor C# in Eq. ©, 
and the parton distributions of Refs. [34,35] with their respective values of strong coupling con- 
stant at LO, NLO and NNLO, a s (M z ) = 0.130, 0.119 and 0.115. Anticipating a slight further 
reduction at N 3 LO, we employ a s (Mz) =0.114 at this order. The N 3 LO corrections to C# in the 
heavy-top limit are taken from Ref. [21]. All higher-order contributions are calculated in the heavy 
top-quark approximation, but normalized to the full lowest-order result. 

As mentioned above, the 8(1 — x) contributions to the N 3 LO coefficient functions c 3 cannot 
be derived at this point. However, we note that the coefficients g Qk in Eqs. (|29l) - (1321 are much 
larger than their 8(1 — x) counterparts for c 1 and especially for c 2 . We expect the same behaviour 
for Cy Moreover, a good approximation (to about 10% or less) to the full double convolutions g* 
g * [ci(x) /x] is obtained at NLO and NNLO by transforming to /V-space and keeping only the \n k N 
and TV terms arising from the +-distributions (but not the 8-function) in c x and c 2 . Consequently 
Eqs. (TT7T ) - (1221) facilitate a sufficient approximation to the complete N 3 LO correction, to which 
we assign a conservative 20% uncertainty, i.e., twice the offset found at NLO and NNLO. 

The corresponding results are added in Fig. [T] to the total cross sections up to NNLO [6-8] at 
the Tevatron and the LHC for the standard choice /j r = ^ = Mh of the renormalization and 
factorization scales. The dependence of the cross sections on yi r is illustrated in Fig. [2] for two 
representative values of the Higgs mass M#. Also shown in Fig. [T] are the additional contribu- 
tions of the N 3 LL threshold resummation (|28T ), see also Ref. [20], of the terms beyond N 3 LO. In 
principle this resummation requires a second coefficient, the four-loop cusp anomalous dimension 
A4, besides D3 of Eq. ([35l ). However, the effect of A4 can safely be expected to be very small, as 
corroborated by the Pade estimate of Ref. [17] employed in our numerical analysis. The Mellin- 
inversion of the exponentiated result (which is entirely dominated by the next few orders in a s in 
the present case) has been performed using the standard 'minimal prescription' contour [36]. 

The inclusion of our new result for the coefficient function c 3 effects an increase of the cross 
sections by about 10% at the Tevatron and 5% at the LHC. The estimated uncertainties due to the 
approximate character of c 3 (see above) thus amount to 2% and 1%, respectively. Contributions of 
yet higher orders, as estimated by the threshold resummation, lead to a further increase by roughly 
half the N 3 LO effect. Lacking N 3 LO and threshold-resummed (see Ref. [37] for a first study) 
parton distributions, the NNLO gluon distribution of Ref. [35] has been employed for all results 
beyond the next-to-leading order. Based on the pattern of the available orders, one may expect 
slightly smaller (by about 2%) gluon-gluon luminosities at N 3 LO. 

The residual uncertainty due to uncalculated contributions of yet higher order is often estimated 
by varying /j r and/or /jjr. At the LHC the representative variation of fj r with fixed /jp illustrated 
in Fig. [2] for two Higgs masses, yields uncertainties of less than 4% for the conventional interval 
1 /2Mh < n r < 2Mh at N 3 LO, an improvement by almost a factor of three over the 9 to 11% sta- 
bility of the NNLO cross sections. At the Tevatron the corresponding /j r dependence decreases 
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100 120 140 160 180 200 100 150 200 250 300 
M H (GeV) M H (GeV) 

Figure 1 : The perturbative expansion of the total cross section for Higgs boson production at the 
Tevatron (left) and the LHC (right) for the standard scale choice fj r =/jjr= M H . 
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Figure 2: The dependence of the fixed-order predictions for the LHC cross section on the renor- 
malization scale /j r at /jj- = Mh for two representative values of the Higgs boson mass M# . 
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from about 11% at NNLO to 5% at N 3 LO for M H = 120 GeV where, as in Fig. El the N 3 LO 
cross section exhibits a stationary point close to /j r = 1/2M#. Considering these and the above 
results, 5% at the LHC and 7% at the Tevatron appear to represent conservative estimates of the 
improved cross-section uncertainties due to the truncation of the perturbation series at N 3 LO. 
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Note added: Shortly after the completion of this letter, Ref. [40] appeared, which addresses the 
threshold resummation especially for lepton-pair production in the approach of Ref. [33]. In par- 
ticular, our result (1331 ) for the coefficient D3 for the Drell-Yan process is confirmed by this research. 
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